Can We Rigorously Define Phases in a Finite System? 



K. A. Bugaev 1 , A. I. Ivanytskyi 1 , E. G. Nikonov 2 , A. S. Sorin 3 and G. M. Zinovjev 1 
1 Bogolyubov Institute for Theoretical Physics, Kiev 03680, Ukraine 
2 Laboratory for Information Technologies, JINR Dubna, 141980 Dubna, Russia 
3 Bogoliubov Laboratory of Theoretical Physics, JINR Dubna, 141980 Dubna, Russia 

(Dated: June 30, 2011) 

Here we propose the generalized statistical multifragmentation model which includes the liquid 
phase pressure of the most general form. This allows us to get rid of the absolute incompressibility of 
the nuclear liquid. Also the present model employes a very general form the surface tension coefficient 
of nuclear fragments. Such a model is solved analytically for finite volumes by the Laplace-Fourier 
transform method to isobaric ensemble. A complete analysis of the isobaric partition singularities 
of this model is also done for finite volumes. It is shown that the real part of any simple pole 
of the isobaric partition defines the free energy of the corresponding state, whereas its imaginary 
part, depending on the sign, defines the inverse decay /formation time of this state. The developed 
formalism allows us to exactly define the finite volume analogs of gaseous, liquid and mixed phases 
of the class of similar models from the first principles of statistical mechanics and demonstrate the 
pitfalls of earlier works. The finite width effects for large nuclear fragments and quark gluon bags 
are also discussed. 

PACS numbers: 25.70. Pq, 21.65.+f, 24.10. Pa 

I. INTRODUCTION 

The necessity to extend the theory of 1-st order liquid-gas phase transition (PT) is determined both by an academic 
interest to this problem and by practical purpose to study the phase transformations in the systems that do not have 
thermodynamic limit. Such an extension is of particular interest for nuclear physics of intermediate energies where 
the nuclear liquid-gas PT [IH3] is studied in the presence of the Coulomb interaction. Also recently there appeared a 
great attention to PTs in finite systems because of the searches for the new state of matter, the quark gluon plasma, 
and its (tri)critical endpoint in the relativistic collisions of heavy ions [2HS]. 

This problem has a long history, but up to now it is not resolved. One of the first attempts of its resolving was 
made by T. Hill [7] whose approach is based on the formulation of thermodynamics of small systems. Hills ideas were 
developed further in [5], where the authors claimed to establish a one-to-one correspondence between the bimodal 
structure of the partition of measurable quantity B, known on average, and the properties of the Lee-Yang zeros [9] 
of this partition in the complex g-plane. However, the analysis of such a definition performed in [10] on the basis of 
exactly solvable model for finite volume [TT] shows that the definition of phases suggested in [5] cannot be established 
experimentally under any circumstances. Therefore, although many important aspects of the nuclear liquid-gas PT 
related to finite volumes of studied nuclear systems are well understood [IJ|3j[l2] the systematic and rigorous extension 
of the PT theory to finite systems is just at the very initial stage in spite of too optimistic beliefs of Ref. [T5] . 

Therefore, here we would like to discuss a powerful mathematical method invented recently [TT] , the Laplace-Fourier 
transform (LFT), that allows one not only to solve analytically several statistical models for finite volumes or surfaces, 
but also to establish a common framework to study the deconfinement PT and nuclear liquid-gas PT in finite systems. 
LFT was successfully applied to the simplified version of the statistical multifragmentation model (SMM) Q31 [TS] in 
finite volumes jllj , to the analysis of surface partition and surface entropy of large but finite physical clusters for a 
variety of statistical ensembles [TS] and to exact solution of the gas of quark gluon bags [17 for finite volumes |10j . 
Furthermore, using the theorems proven in |10[ 118] it is possible to straightforwardly apply the exact representation 
of the finite volume grand canonical partition (GCP) of the gas of quark gluon bags model [TU] to other exactly 
statistical models of the deconfinement PT which were solved recently in thermodynamic limit [19H23] . Finally, we 
would like to stress that our approach to study PTs in finite systems is not restricted to the models without long- 
range Coulomb-like interaction. Our strategy is as follows: to firmly define the phases in finite systems without the 
long-range interaction, and then to extend this approach to the systems with the Coulomb-like interaction. This work 
is mainly devoted to the first of these tasks for which the finite volume solution of the simplified version of the SMM 
(and similar models) is analyzed here. In addition, here we generalize the simplified version of the SMM [T4J [15] in 
order to repair such its defects as an existence of limiting baryonic density and an absence of (tri) critical endpoint for 
the values of the Fisher exponent r > 2 [2~4"l 125] . 

The work is organized as follows. Section II is devoted to the formulation of the generalized SMM. Section III 
contains an introduction into the LFT technic. The analysis of the singularities of the isobaric partition of the 
suggested model is given in Section IV. In Sections V and VI we, respectively, discuss the location of the isobaric 
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partition singularities in complex plane for the case without PT and with PT in thermodynamic limit. The finite 
volume analogs of phases along with the critique of the Hill's treatment of 1-st order PT in finite systems are presented 
in Section VI. Finally, Section VII contain our concluding remarks. 



II. GENERALIZED SMM 



The system states in the standard SMM are specified by the multiplicity sets {n^} {n^ = 0, 1, 2, ...) of fc-nucleon 

fragments. The partition function of a single fragment with k nucleons is [1]: V(j)k{T) = V (rnTk/2n) 3 ^ 2 Zk , where 
k = 1,2,...,A (A is the total number of nucleons in the system), V and T are, respectively, the volume and the 
temperature of the system, m is the nucleon mass. The first two factors on the right hand side (r.h.s.) of the single 
fragment partition originate from the non-relativistic thermal motion and the last factor, Zk, represents the intrinsic 
partition function of the k- nucleon fragment. Therefore, the function <pk{T) is a phase space density of the k- nucleon 
fragment. For k = 1 (nucleon) we take Z\ = 4 (4 internal spin-isospin states) and for fragments with k > 1 we use 
the expression motivated by the liquid drop model (see details in Ref. [T]): Zk = exp(— fk/T), with the fragment free 
energy of the simplified SMM 

fk = -WIT) k + o(T) 3 + (t + 3/2)T In k , (1) 

with W(T) = W a + T 2 /e Q . Here W a — 16 MeV is the bulk binding energy per nucleon. T 2 /e D is the contribution of 
the excited states taken in the Fermi-gas approximation (e = 16 MeV). er(T) is the temperature dependent surface 
tension parameterized in the following relation: o(T) = a a [{T 2 - T 2 )/(T 2 + T 2 )] 5 / 4 , with a a = 18 MeV and 
T c = 18 MeV (ct = at T > T c ). The last contribution in Eq. (JlJ involves the Fisher's term with dimensionless 
parameter r. The free energy ([!]) does not contain the symmetry and Coulomb contributions which are neglected. 
Hence it is called a simplified version of the SMM which was suggested and studied numerically in Refs. [TH I15j. 
However, its investigation appears to be of principal importance for studies of the liquid-gas phase transition in finite 
systems. 

From an exact analytical solution of the simplified SMM [53] it follows that the baryonic density has a maximal 
value max/) = 1/6 = p Q = 0.16 fm -3 in the limit [i — > oo (here p Q is the normal nuclear density). This, however, 
contradicts to the experiments on heavy ion collisions |26j in which the nuclei can be compressed to much higher 
densities. In order to avoid such a behavior we generalize the SMM to the GSMM which includes the pressure of the 
liquid phase pi(T,fi) in a general form. Then the free energy of the fragments with k > 1 reads as 

fjf = pk-pi(T,fj,)bk + a(T) k 2/3 + (r + 3/2)Tln k , (2) 

where the pressure of a liquid phase is (at least) a double diffcrcntiable function of its arguments that contains the 
temperature dependent binding energy W (T). Also it is assumed that the function pi(T,p) reproduces all the typical 
properties of a liquid phase. An extremely important property of the GSMM is that the liquid phase equation of state 
Pi(T,p) can be taken from some microscopic models including the mean-field ones, but the resulting model will be 
a truly statistical one since the analytical properties of the isobaric partition singularities remain unmodified in this 
case. 

The GSMM nucleons are considered as in the SMM. Note that the pressure of a liquid phase should approach 
the asymptotics pi ~ T 2 for T — > oo and pi ~ fj, 2 for /i —> oo to respect a causality condition [27]. Therefore, the 
simplest parameterization of the liquid phase pressure which at low densities recovers the usual SMM result and at 
high densities obeys such asymptotics can be written as follows 

[i(l + ay) + W(T) 

Pl{T,p,)= 7 , (3) 

where a positive constant a > has to be fixed by the condition that at low densities it behaves as a\/i\ <C 1. However, 
in what follows we study the most general form of the liquid phase pressure. 

In addition to the new parameterization of the free energy of the fc-nucleon fragment ^ we propose to consider a 
more general parameterization of the surface tension coefficient 



a(T) 



T r -T 



T, 



sign(T c - T) , (4) 



with C = const > 1 and T c = 18 MeV. In contrast to the Fisher droplet model [28] and the SMM [T], the GSMM 
surface tension Q is negative above the critical temperature T c . It is necessary to stress that there is nothing wrong 
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or unphysical with the negative values of surface tension coefficient (|4|), since a k 3 in (|2|) is the surface free energy of 
the fragment of mean volume b k and, hence, as any free energy, it contains the energy part e sur f and the entropy 
part s sur f multiplied by temperature T [28 . Therefore, at low temperatures the energy part dominates and the 
surface free energy is positive, whereas at high temperatures the number of fragment configurations with large surface 
drastically increases and it exceeds the Boltzmann suppression and, hence, the surface free energy becomes negative 
since s sur f > e °y / . Because of this reason the negative values of the surface tension coefficient were recently employed 
in a variety of exactly solvable statistical models for the deconfinement PT [L9Tf21~1 [29] . For the first time this fact 
was derived within the exactly solvable models for surface deformations of large physical clusters [16] . Very recently 
two of us derived a relation between the surface tension of large quark gluon bags and the string tension of two static 
color charges measured by the lattice QCD [30] from which it was possible to conclude that at high temperatures the 
surface tension coefficient of quark gluon bags should be negative [30, 31,. 

Furthermore, a thorough analysis of the temperature dependence of the surface tension coefficient in ordinary 
liquids [3U [33J shows not only that the surface tension coefficient approaches zero, but, in contrast to the widely 
spread beliefs, for many liquids the full T derivative of a does not vanish and remains finite at T c : ff < [3j. 
Therefore, just the naive extension of these data to the temperatures above T c would lead to negative values of surface 
tension coefficient at the supercritical temperatures. On the other hand, if one, as usually, believes that a = for 
T > T c , then it is absolutely unclear what physical process can lead to simultaneous existence of the discontinuity 
of 4-^ at T c and the smooth behavior of the pressure's first and second derivatives at the cross-over. Therefore, we 
conclude that negative values of the surface tension coefficient at supercritical temperatures are also necessary for 
ordinary liquids although up to now this question has not been investigated. 

The existence of negative values of the surface tension coefficient in ^ leads to entirely new result for the GSMM 
compared to that one of the SMM for t > 2. Thus, for r > 2 the SMM predicts an existence of the 1-st order PT 
up to infinite values of T [TH] [231 Hi] . Clearly, such a result does not correspond to the experimental findings and is 
usually understood as a pitfall of this model. However, the negative values of a in Q lead to a different result in the 
GSMM. Using the technic developed in [13 [5U] it is easy to show that in this case there is a cross-over for T > T c 
and, hence, for r > 2 the GSMM has a critical point at T = T C . 



III. THE LAPLACE-FOURIER TRANSFORMATION TECHNIC 



To evaluate the GCP of the GSMM for finite volumes first we define the canonical partition function (CPF) of 
nuclear fragments. The latter has the following form: 



Zf{V,T) 



E 



n 

fc=i 



[y MT)y 



(5) 



In Eq. ([5]) the nuclear fragments are treated as point-like objects. However, these fragments have non-zero proper 
volumes and they should not overlap in the coordinate space. In the excluded volume (Van der Waals) approximation 
this is achieved by substituting the total volume V in Eq. 151 by the free (available) volume Vf = V — bJ2k knt, 
where b is eigen volume of nucleon. Therefore, the corrected CPF becomes: Za(V,T) = Z l /(V — bA,T). 

The calculation of Za{V, T) is difficult due to the constraint J2k = A. This difficulty can be partly avoided by 
evaluating the GCP function: 



00 

Z(V,T,a) = ]Texp(^)z A (V,T) Q(V - bA) 



(6) 



A=0 



where n denotes a chemical potential. Nevertheless, the calculation of Z is still rather difficult. The summation over 
{nfe} sets in Za cannot be performed analytically because of additional A-dependence in the free volume Vf and the 
restriction Vf > 0. This problem was resolved [23] by the Laplace transformation method to the so-called isobaric 
ensemble [17] . 

K(V) 

To study the PT in finite systems here we consider a more strict constraint ^ k = A , where the size of the 

k 

largest fragment K (V) = aV/b cannot exceed the total volume of the system (the parameter a < 1 is introduced for 
convenience). The case K(V) = const considered in [34] is also included in our treatment. A similar restriction should 
be also applied to the upper limit of the product in all partitions Z^(V, T), Za{V, T) and Z(V, T, u) introduced above 
(how to deal with the real values of K(V), see [H]). Then the model with such a constraint, the CGSMM, cannot 
be solved by the Laplace transform method, because the volume integrals cannot be evaluated due to a complicated 
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functional ^-dependence. However, the CGSMM can be solved analytically with the help of the following identity 

eld 



Hoo +00 



G{V) 



dr\ 
2^ 



Mv-6) 



G(0, 



(7) 



—00 —00 



which is based on the Fourier representation of the Dirac <5-function. The representation ^ allows us to decouple 
the additional volume dependence and reduce it to the exponential one, which can be dealt by the usual Laplace 
transformation in the following sequence of steps 



+00 +00 



Z(X,T,fi) 



dV e- xv Z{V,T,n) = J dV' J d£ 

— OO — OO 

ff £ {V MT)e-=^} 

k—1 

-\-oc -\-oc 

I?' hit 



p iv{V'-£.)-W' 

2tt 



E 

{rik} 



e(V) 



Jn(V'-S)-\V'+V'F(t,X-iri) 



(8) 



-OO — OO 



After changing the integration variable V ^ V = V — b k n k, the constraint of 0-function has disappeared. Then 

k 

all rtfc were summed independently leading to the exponential function. Now the integration over V in Eq. (|) can 
be straightforwardly done resulting in 



hoo +00 



Z{X,T,n)= dt; 



dij 

2ir X — if] — J-(^, X — irj) 



-00 —00 



where the function X) is defined as follows 



K{0 
k=l 



-\bT)k ( mT\ 



K(0 

(i XbT \ -* _ 

zi e t +2_^ k t 

k=2 



{ Pl (T,ii)-\T)bk-<7k* 



(9) 



(10) 



This result generalizes the finite volume solution of the simplified SMM obtained in [TUl QT] . 

As usual, in order to find the GCP by the inverse Laplace transformation, it is necessary to study the structure of 
singularities of the isobaric partition (10). 



IV. ISOBARIC PARTITION SINGULARITIES 



The isobaric partition ( 10 1 of the CGSMM is, of course, more complicated than its SMM analog found in thermo- 



dynamic limit |24] because for finite volumes the structure of singularities in the CGSMM is much richer than in the 
SMM, and they match in the limit V — > 00 only. To see this let us first make the inverse Laplace transform: 



X+ioo 



Z(V,T,n) 



^.Z{X,T,») 



e \V-i V (. 



+00 +00 X+ io ° 

r dn r dA 

4 J 2?r J 2ni X - ir\ - J"(£> X - irj) 

— oo x~i°° 



3J r (g,A„) 
3A„ 



(11) 



where the contour A-integral is reduced to the sum over the residues of all singular points X = X n + ir] with n = 1, 2, .., 
since this contour in the complex A-plane obeys the inequality \ > max(i?e{A n }). Now both remaining integrations 
in (11) can be done, and the GCP becomes 



Z(V,T,n) 



E e 



X n V 



1 



(12) 
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i.e. the double integral in (11) simply reduces to the substitution £ — > V in the sum over singularities. In [11] this 



remarkable result was formulated as a theorem, which now is generalized to more complicated forms of the liquid 
phase pressure. 



The simple poles in (11 1 are defined by the equation 



A„ = JYV,A n ) 



(13) 



In contrast to the usual SMM [21] the singularities A„ are (i) are volume dependent functions, if K(V) is not constant, 
and (ii) they can have a non-zero imaginary part, but in this case there exist pairs of complex conjugate roots of ( 13 ) 
because the GCP is real. 

Introducing the real R n and imaginary /„ parts of A n = R n + il n , we can rewrite Eq. (13 1 as a system of coupled 
transcendental equations 



K(V) 

k=l 
K(V) _ 

I n = - ^ 4> k {T) e ' r sm(I n bk) , 

k=l 

where we have introduced the set of the effective chemical potentials v n = v(\ n ) with v(X) 



\ 2tt / 



* Ziexp((jJt-pi(T,iJt)b)/T) and 4> k >i(T) = ( 
0), first. Similarly to the SMM, for /„ = / 



2-rt ) 



(14) 
(15) 

Pi(T, /i)6 — A6T, and 
T exp(-a(T) k 2 / 3 /T) 



the real root i? n of the 



the reduced distributions 4>i(T) 
for convenience. 

Consider the real root (Rq > 0, Iq 
GSMM exists for any T and fi. Comparison of i?o from ( 14 1 with the expression for vapor pressure of the analytical 
SMM solution [21] indicates that TRq is a constrained grand canonical pressure of the mixture of ideal gases with 
the chemical potential vq. As usual, for finite volumes the total mechanical pressure [Zj differs from TRq. Equation 
( [15] ) shows that for I n> o ^ the inequality cos(I n bk) < 1 cannot simultaneously become an equality for all fc-values. 
Then from Eq. ( 14 1 one obtains (n > 0) 



R, 



K(V) 
k=l 



Rn < Ro ■ 



(16) 



where the second inequality ( 16 ) immediately follows from the first one. In other words, the gas singularity is always 
the rightmost one. This fact plays a decisive role in the thermodynamic limit V — > 00. 

The interpretation of the complex roots A rl> o seems to be less straightforward and, hence, in this case we follow 
the line of arguments suggested in Ref. [Hi. According to Eq. (12), the GCP is a superposition of the states of 
different free energies — X n VT. Strictly speaking, — X n VT has a meaning of the change of free energy, but we will use 
the traditional term for it. For n > the free energies are complex. Therefore, —\ n >oT is the density of free energy. 
The real part of the free energy density, —R n T, defines the significance of the state's contribution to the partition: 



due to ( 16 ) the largest contribution always comes from the gaseous state and has the lowest value of the real part of 



free energy density. As usual, the states which do not correspond to the lowest value of the (real part of) free energy, 
i. e. —R n> oT, are thermodynamically metastable. For infinite volume they should not contribute unless they are 
infinitesimally close to —RqT, but for finite volumes their contribution to the GCP may be important. 



As one can see from ( 14 ) and ( 15 ), the states of different free energies have different values of the effective chemical 
potential v n , which is not the case for infinite volume [H], where there exists a single value for the effective chemical 



potential. Thus, for finite V the states which contribute to the GCP (12) are not in a true chemical equilibrium 



The meaning of the imaginary part of the free energy density becomes clear from ( 14 ) and ( 15 1 : as it is seen from ( 14 ) 



the imaginary part I rl >o effectively changes the number of degrees of freedom of each fc-nucleon fragment (k < K{V)) 
contribution to the free energy density ~R n> oT. It is clear, that the change of the effective number of degrees of 
freedom can occur virtually only and, if A„>o state is accompanied by some kind of equilibration process. Both of 
these statements become clear, if we recall that the statistical operator in statistical mechanics and the quantum 
mechanical evolution operator are related by the Wick rotation [35] . In other words, the inverse temperature can be 
considered as an imaginary time. Therefore, depending on the sign, the quantity I n bT = r^ 1 that appears in the 



trigonometric functions of the equations (14) and (15) in front of the imaginary time 1/T can be regarded as the 
inverse decay/formation time r„ of the metastable state which corresponds to the pole A„>o (for more details see next 
sections and [TT]). 

Such an interpretation of t„ naturally explains the thermodynamic metastability of all states except the gaseous 
one: the metastable states can exist in the system only virtually because of their finite decay /formation time, whereas 
the gaseous state is stable because it has an infinite decay /formation time. 
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FIG. 1: A graphical solution of Eq. (15 1 for T — 10 MeV and r = 1.825 for the typical SMM parameterization of the surface 
tension coefficient by Eq. |TJ. Note, however, that qualitatively the same picture remains valid for any parameterization of 
the surface tension coefficient. The l.h.s. (straight line) and r.h.s. of Eq. ( 15 1 (all dashed curves) are shown as the function of 
dimensionless parameter Ii b for the three values of the largest fragment size K(V). The intersection point at (0; 0) corresponds 
to a real root of Eq. 1 13). Each tangent point with the straight line generates two complex roots of 1 13 1. 



V. NO PHASE TRANSITION CASE 



It is instructive to treat the effective chemical potential i^(A) as an independent variable instead of fi. In contrast to 
the infinite V, where the upper limit v < defines the liquid phase singularity of the isobaric partition and gives the 
pressure of a liquid phase pi(T, fx) = TRq\v^oo, for finite volumes and finite K(V) the effective chemical potential can 
be complex (with either sign for its real part) and its value defines the number and position of the imaginary roots 
{A„>o} in the complex plane. Positive and negative values of the effective chemical potential for finite systems were 
considered within the Fisher droplet model |36j . but, to our best knowledge, its complex values have been discussed 
for the first time in [TT]. From the definition of the effective chemical potential v(X) it is evident that its complex 
values for finite systems exist only because of the excluded volume interaction, which is not taken into account in the 
Fisher droplet model [28] . 

As it is seen from Fig. 1, the r.h.s. of Eq. (15) is the amplitude and frequency modulated sine-like function of 
dimensionless parameter /„ b. Therefore, depending on T and Re(v) values, there may exist no complex roots {A n >o}, 
a finite number of them, or an infinite number of them. In Fig. 1 we showed a special case which corresponds to 
exactly three roots of Eq. ( 13 1 for each value of K(V): the real root (Iq — 0) and two complex conjugate roots (±-Tl). 
Since the r.h.s. of (15) is monotonously increasing function of Re{v) 1 it is possible to map the T — Re(v) plane into 
regions of a fixed number of roots of Eq. (13). For fixed T- value each curve in Fig. 2 divides the T — Re(v) plane 
into three parts: for i?e(^)-values below the curve Re(vi(T)) there is only one real root (gaseous phase), for points 
on the curve Re(v) = Re{v\{T)) there exist three roots, and above the curve Re(vi(T)) there are five or more roots 
of Eq. (13). Although Fig. 2 corresponds to the usual SMM parameterization of the surface tension coefficient, the 
picture is qualitatively the same for general parameterization of a whereas its modifications are discussed below. 

For constant values of K(V) = K the number of terms in the r.h.s. of ( |15[ ) does not depend on the volume and, 
consequently, in thermodynamic limit V — > oo only the rightmost simple pole in the complex A-plane survives out 
of a finite number of simple poles. According to the inequality (16), the real root Ao is the rightmost singularity of 
isobaric partition (J9J). However, there is a possibility that the real parts of other roots A n> o become infinitesimally 
close to Rq, when there is an infinite number of terms which contribute to the GCP (12). 

Let us show now that even for an infinite number of simple poles in ( 12 ) only the real root Ao survives in the limit 
V —> oo. For this purpose consider the limit Re{v n ) 3> T. In this limit the distance between the imaginary parts of 
the nearest roots remains finite even for infinite volume. Indeed, for Re(y n ) 3> T the leading contribution to the r.h.s. 
of ( 15 ) corresponds to the harmonic with k — K, and, consequently, an exponentially large amplitude of this term 
can be only compensated by a vanishing value of sin (/„ bK), i.e. /„ bK = im + 6 n with \5 n \ <C vr (hereafter we will 
analyze only the branch /„ > 0), and, therefore, the corresponding decay/formation time T n « K^nT]" 1 is volume 
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FIG. 2: Each curve separates the T — Re(v n ) region of one real root of Eq. ( |13| l (below the curve), three complex roots (at the 
curve) and five and more roots (above the curve) for three values of K (V) and the same parameters as in Fig. 1. 



independent. 



Keeping the leading term on the r.h.s. of (15) and solving for 8 n , one finds 

7T n + 5 n 



In 



bK{v) ' 

(-l)" + 1 7m gejaQK 

Rn « (-l) B ^Jc(r) e 1 ^^ , 



where in the last step we used Eq. (14) and condition \S n \ <C 7r. Since for y — >• oo all negative values of R n 
contribute to the GCP dl2[) , it is sufficient to analyze even values of n which, according to (19), generate R n 

can not be broken, a single possibility, when A„>o 
for some finite n 



Since the inequality 



(17) 

(18) 

(19) 

cannot 
> 0. 



( 12 ), corresponds to the case R n — > R a — (r 



pole can contribute to the partition 
Assuming this, we find Re(v(X n )) — > Re(v(\o)) for the 



same value of \x. 



Substituting these results into equation ( 14 ) , one gets 



R n 



K 

4>k\T) e t cos 

k=i 



nnk 



K 



< Rq. 



(20) 



The inequality (20) follows from the equation for Rq and the fact that, even for equal leading terms in the sums above 

the difference between Rq and R n is large due to the next to leading term k = K — 1, 

He(i.(t )) (K-l) , ... . i 

t 1. Thus, we arrive at a contradiction with our assumption Rq — R n — > + , 



(with k = K and even n) 
which is proportional to e" 



and, consequently, it cannot be true. Therefore, for large volumes the real root Ao always gives the main contribution 
to the GCP ( [12"] ), and this is the only root that survives in the limit V — > 00. Thus, we showed that the model 
with the fixed size of the largest fragment has no phase transition because there is a single singularity of the isobaric 
partition ([9]), which exists in thermodynamic limit. However, for the finite systems we can also define the analog 
of the metastable mixed phase which corresponds to a finite number of complex conjugate solutions A ra> o. Clearly, 
that in thermodynamic limit the contribution of these metastable states into all physical quantities disappears. The 
equation Re(v) = Re{v\{T)) defines the boundary between the finite volume analogs of the gaseous and mixed phases 
in T — Re{v) and T — fi planes. 



VI. FINITE VOLUME ANALOGS OF PHASES 



If KiV) monotonically grows with the volume, the situation is different. In this case for positive value of Re(v) 3> T 
the leading exponent in the r.h.s. of (15) also corresponds to a largest fragment, i.e. to k — K(V). Therefore, we 
can apply the same arguments which were used above for the case K (V) = K = const and derive similarly equations 
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( 17)-(fl9|) for /„ and R n . From the relation /„ 



bK(V) 



it follows that, when V increases, the number of simple poles 



in ( 1 1[) also increases and the imaginary part of the closest to the real A-axis poles becomes very small, i.e /„ — » 
for n <C Kiy), and, consequently, the associated decay /formation time r„ ps R~(V)[irnT]~ 1 grows with the volume 
of the system. Due to the fact that /„ — > 0, the inequality (20) cannot be established for the poles with n <C K(V). 



Therefore, in contrast to the previous case, for large K(V) the simple poles with n <C K(V) will be infinitesimally 
close to the real axis of the complex A-plane. 
In this case from Eq. ( 19 ) one obtains 



TR n ps pi(T,fx) 



Pl(T, fx) 



K{V)b 



In 



[K(V)]h 



<Pk(T) 
-T 



t\tlK{V) 



K{V)b 



(21) 



for Re(v) ^> T and K(V) ^> 1. Thus, from Eq. (21 ) one can clearly see that for an infinite volume an infinite number 
of simple poles moves toward the real A-axis to the vicinity of liquid phase singularity A; = pi(T, fi)/T of the isobaric 
partition |24) and generates an essential singularity of function T{V,pi/T) in (10) irrespective to the sign of the liquid 
phase pressure pi(T, fx) . As we showed above, the states with Re(y) ^> T become stable because they acquire infinitely 
large decay /formation time r n in the limit V — > oo. Therefore, these states should be identified as a liquid phase 
for finite volumes as well. Such a conclusion can be easily understood, if we recall that the partial pressure TR n of 
(21) corresponds to a single fragment of the largest possible size. Moreover, as one can see from the leading terms 
on the r.h.s. of (21) the partial pressure TR n contains both the liquid phase and the surface contributions for a 
spherical fragment of the mean radius [K(V)b}3. In fact, the above results remain valid under a weaker condition 
Re(v)K(V) 3> T since such an inequality allows one to establish the approximation (19). 

Now it is clear that each curve in Fig. 2 is the finite volume analog of the phase boundary T — v for a given value 
of K(V): below the phase boundary there exists a gaseous phase, but at and above each curve there are states which 
can be identified with a finite volume analog of the mixed phase, and, finally, if Re(y)jT — > oo there exists a liquid 
phase. Again as in the previous section the equation Re(v) — Re(v\(T)) defines the boundary between the finite 
volume analogs of the gaseous and mixed phases in T — Re(v) and T — fx planes. Clearly, for finite V the solution of 
this equation fx c (T, V) depends on T and V . 

Although the calculations depicted in Fig. 2 were made for a(T) > and for finite values of the effective chemical 
potential Re(v) the shown results can be qualitatively explained using Eq. (21) in the limit Re(v)K(V)/T — > oo. 
Indeed, from |2ll) one finds 



Re{u n ) 



[K(V)]i 



T 



ln|( 



2tt *i 



R n \+TlnK{V) 



K(V) 



(22) 



from which one obtains Re{v n ) > for finite K(V) values (compare to Fig. 2), if | (^^) 2 R n \ 3> 1 and a > 0. If, 
however, the surface tension coefficient gets negative, i.e. a < for T > T c , then for sufficiently large values of K(V) 
one can find that Re{v n ) < 0, i.e. in this case the finite volume analog of phase boundary can demonstrate another 
behavior than that one shown in Fig. 2. For T = T C the surface tension coefficient vanishes and from (22) we get 



Re{v n ) 



T=T C 



K(V) 



In 



2tt y pi(T c ,tx) 
rnTr J T r 



■T\nK(V) 



(23) 



where in the last step of derivation we replaced R n by the leading term from the r.h.s. of (21). This result shows that 



(i) at T = T c the deviation of the partial pressure TR n from the liquid phase pressure decreases faster as function 
of Kiy) than for other temperatures, but at the same time (ii) for large K(V) this deviation still decreases slower 
than the imaginary part /„. For a quantitative example let us choose /x in such a way that the first term on the r.h.s. 
of (23) disappears, i.e. for Re(v n )\ T=T « -^f^h\K{V). Then for r = 1.825 one finds Re(v n )\ T=T ~ 7.56 MeV for 

3.72 MeV for K{V) = 30, and Re{v n )\ T=T ^ 



It=t c ~ k{v) 

K{V) = 10, Re{v n )\ T=T ^ 4.92 MeV for K{V) = 20, Re(u n )\ T=T( 
MeV for KCV) = 100. From Fig. 2 one can see that, although our estimate of Re{u n )\ for KCV) — 10 is about 
2.5 MeV below its value found numerically, the corresponding estimates for K(V) — 20 and K(V) — 30 obtained from 
(23) are in a very good agreement with the results of the numerical evaluation. 

When there is no phase transition, i.e. K (V) = K = const, the structure of simple poles is similar, but, first, the 
line which separates the gaseous states from the metastable states does not change with the system volume, and, 
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second, as shown above, the metastable states will never become stable. Therefore, a systematic study of the volume 
dependence of free energy (or pressure for very large V) along with the formation and decay times may be of a crucial 
importance for experimental studies of the nuclear liquid gas phase transition. 

The above results demonstrate that, in contrast to Hill's expectations [7J, the finite volume analog of the mixed 
phase does not consist just of two pure phases. The mixed phase for finite volumes consists of a stable gaseous phase 
and the set of metastable states which have different free energy values. Moreover, the difference between the free 
energies of these states is not the surface- like, as Hill assumed in his treatment [7J, but the volume-like as we have 



seen. Furthermore, according to Eqs. (14) and (15 1, each of these states consists of the same fragments, but with 
different weights. As was shown above for the case Re(v) ^> T, some fragments that belong to the states, in which 
the largest fragment is dominant, may, in principle, have negative weights (effective number of degrees of freedom) 



in the expression for R n >o (14). This can be understood easily because higher concentrations of large fragments can 
be achieved at the expense of the smaller fragments and is reflected in the corresponding change of the real part of 
the free energy ~R u> qVT. Therefore, the actual structure of the mixed phase at finite volumes is more complicated 
than was expected in earlier works. 

Here it is necessary to add a few remarks about the description of the deconfinement PT on the basis of statistical 
models 125] that were solved recently in thermodynamic limit. The finite volume solution of the models [T51 - 

[2T] [25J can be straightforwardly found using the LFT developed in jTHl [TT] . Recently, however, the importance of 
finite width of heavy/large quark gluon plasma bags was realized [2TM23] . Both the theoretical estimates [22] [23] and 
the analysis of the asymptotic Regge trajectories of non-strange mesons [37] indicate that the width of the quark gluon 

bag of the volume V with the mass M being heavier than M w 2.5 GeV is T = 7(T)y^ (here V > V Q = 1 fm 3 ). 

Since even at T = the value j(T = 0) « 400 MeV is large, we conclude that the short life time tuf e (V) = l/r(V) of 
such bags can, in principle, modify our conclusions about the metastable states A„>o- This is so because in the SMM 
and the GSMM the fragments are implicitly assumed to be stable and, hence, their life-time is set to be infinite while 
in reality there is the region of stability of nuclear fragments outside of which the nucleus life-time is extremely short. 
A similar situation is with large/heavy quark gluon bags although the instability of nuclei is due to the Coulomb 
interaction whereas the short life-time of the bags should be attributed to such a property of strongly interacting 
matter as the color confinement [2T1 [3U1 13T] . 

Now it is clear that, if the individual life-time of the largest bag tuf e (V) in the volume V is larger than the collective 
decay/formation time of the state A n >o which for quark gluon bags is r„ f=a Vg % nT [10] . then the largest bag can be 
considered as a stable one during the course of collective decay /formation process, i.e. tuf e (V) > r n . Otherwise the 
process of collective decay or formation of the state A„>o cannot ever be completed because the largest bag ceases its 
existence much earlier. Note that under such conditions one can hardly expect an existence of thermal equilibrium 
in the system of shortly living particles. 

Therefore, the inequality tuf e (V) > t„ sets some constraints on the applicability range of the above analysis to the 
short-living bags 



< ^Ln. (24) 



'V 

This inequality can be used as an estimate for the volume of largest bags in a finite system, if n a nd T are known 



For instance, if 7(T) = j(T = 0) = const, then one obtains X- < 



7 (T=0) 



Alternatively, Eq. (24) can be used 



to estimate the range of temperatures for given n since V/Vq > 1. Using the last inequality one can rewrite (24) as 
n > Substituting into the last result the following expression for j(T) = ^2V Q aTTj 1 (T 2 + TT H + Tf) that 

was predicted in [22], where T# m 160 MeV is the Hagedorn temperature and a = |^7r 2 for SU(3) color group with 
two quark flavors, we find n > 1 for T > 500 MeV and n > 2 for T > 1120 MeV. These estimates show that at LHC 
energies we can expect an existence of metastable states which can be described within the developed approach. 



VII. CONCLUSIONS 



In this work we generalized the SMM to GSMM and included into it the nuclear liquid phase pressure of the most 
general form. This allowed us to get rid of the absolute incompressibility of the nuclear liquid. Also in this model 
we introduced very general form of the surface tension coefficient which enabled us to repair another pitfall of the 
simplified SMM related to the absence of critical endpoint in this model for r > 2. The LFT method was applied to 
the constrained GSMM and this model was solved analytically for finite volumes. 

It is shown that for finite volumes the GCP function can be identically rewritten in terms of the simple poles of 
the isobaric partition fcty . The real pole Aq exists always and the quantity TXq is the GCP pressure of the gaseous 



10 



phase. The complex roots A ra >o appear as pairs of complex conjugate solutions of equation (13). As we discussed, 
their most straightforward interpretation is as follows: — T Re(X n> o) has a meaning of the free energy density, whereas 
bTIm(X n> o), depending on sign, gives the inverse decay /formation time of such a state. The gaseous state is always 
stable because its decay/formation time is infinite and because its free energy is minimal. The complex poles describe 
the metastable states for i?e(A„>o) > and mechanically unstable states for i?e(A n>0 ) < 0. 

We studied the volume dependence of the simple poles and found a dramatic difference in their behavior in case 
with PT and without it. For the former one the found representation allows one to define the finite volume analogs 
of phases unambiguously and to establish the finite volume analog of the T — fi phase diagram (see Fig. 2). At 
finite volumes the gaseous phase exists, if there is a single simple pole, the mixed phase corresponds to three and 
more simple poles, whereas the liquid is represented by an infinite amount of simple poles at highest possible particle 
density (or /i — > oo). 

As we showed, for given T and fi the states of the mixed phase which have different i?e(A„) are not in a true chemical 
equilibrium for finite volumes. This feature cannot be obtained within the Fisher droplet model due to lack of the 
hard core repulsion between fragments. This fact also demonstrates clearly that, in contrast to Hill's expectations 
[7], the mixed phase is not just a composition of two states which are the pure phases. As we showed, the mixed 
phase is a superposition of three and more collective states, and each of them is characterized by its own value of A n . 
Because of that the difference between the free energies of these states is not the surface- like, as Hill argued [7j, but 
the volume-like. 

A detailed analysis of the isobaric partition singularities in the T — Re(v) plane allowed us to define the finite 
volume analogs of phases and study the behavior of these singularities in the limit V — > oo. Such an analysis opens 
a possibility to rigorously study the nuclear liquid-gas PT and the deconfinement PT directly from the finite volume 
partition. This may help to extract the phase diagram of strongly interacting matter from the experiments on finite 
systems (nuclei) with more confidence. The conditions of the model applicability to the description of the short-living 
quark gluon plasma bags are also discussed. 
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